Polarized fermions in the unitarity limit 
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We consider a polarized Fermi gas in the unitarity limit. Results are calculated analytically up to next-to- 
leading order in an expansion about d = 4 spatial dimensions. We find a first order transition from superfluid 
to normal phase. The critical chemical potential asymmetry for this phase transition is 6/./c = ^(1 — 0.467e), 
where e = 4 — is the expansion parameter and /j is the average chemical potential of the two fermion species. 
Stability of the superfluid phase in the presence of supercurrents is also studied. 
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I. INTRODUCTION 



Recently, there has been a lot of interest in the quantum 
phase transition between a paired fermion superfluid and a 
normal Fermi liquid that occurs as the difference in the chem- 
ical potentials of the up and down spins is increased. This 
transition is well understood in the strong coupling Bose- 
Einstein condensation (BEC) and weak coupling Bardeen- 
Cooper-Schrieffer (BCS) limits, but the nature of the phase di- 
agram near the BEC/BCS crossover remains to be elucidated. 
The BEC/BCS crossover is characterized by a divergent atom- 
atom scattering length. Over the last year, first results from 
experiments with polarized atomic systems near a Feshbach 
resonance have appeared Ol 0.1111 ■ 

Theoretically the regime of large scattering lengths is dif- 
ficult since the standard perturbative methods are not appli- 
cable. Based on an observation by Nussinov and Nussinov 
101, Nishida and Son fsl proposed an analytic method for cal- 
culating thermodynamic properties as an expansion around 
d = 4- spatial dimensions. One starts by performing the cal- 
culation in arbitrary li = 4 — e space dimensions and develops 
a perturbative expansion in e. In this formalism, the Cooper 
pair energy Xq is considered 0(1) and the chemical potential 
/7 ^ 0(e). A next-to-leading order calculation jsl of the su- 
perfluid gap and the equation of state agrees well with fixed 
node Green Function Monte Carlo |0, US 1^ and Euclidean 
Path Integral iflolfTTIl calculations. Somewhat smaller values 
of the energy per particle and the energy gap were obtained 
in the canonical Path Integral Monte Carlo calculation lfl2ll . 
For a polarized system with chemical potential difference 5/7, 
we expect the superfluid phase to become unstable when the 
asymmetry is on the order of the gap in the symmetric system 
Sa* - Xo ill HlISI. Thus we will consider a situation where 
A* ~ e < Xo ~ 1 ~ 5/7- 
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II. EPSILON EXPANSION 

The physics of the unitarity limit is captured by an effective 
lagrangian of point-like fermions interacting via a short-range 
interaction. The lagrangian is 



L = \\i 




(1) 



where i|/ = (i|/|,\|/j^) is a two-component spinor. The coupling 
constant Co is related to the scattering length. In dimensional 
regularization the unitarity limit a ^ °o corresponds to Co 
oo. In this limit the fermion-fermion scattering amplitude is 
given by 
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where 8 0+. As a function of d the Gamma function has 
poles at d — 2,4, . . . and the scattering amplitude vanishes at 
these points. Near d ^ 2 the scattering amplitude is energy 
and momentum independent. For d = 4 — ewe find 



STt^e 



PQ — Ep/l + ih 



(3) 



We observe that at leading order in e the scattering ampli- 
tude looks like the propagator of a boson with mass 2m. The 
boson-fermion coupling is ~ (Sti's) / rn^ and vanishes as 
e ^ 0. This suggests that we can set up a perturbative expan- 
sion involving fermions of mass m weakly coupled to bosons 
of mass 2m. We can ehminate the four-fermion coupling in 
Eq. ([T) using a a Hubbard-Stratonovich transformation. In the 
unitarity limit we get 
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ido + h^i + Oi — 
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^^a+1'(j) + /i.c. 



(4) 



where *P = (\|/|,\|/|)^ is a two-component Nambu-Gorkov 
field, a,- are Pauli matrices acting in the Nambu-Gorkov space 
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and a± = (ai ± i02) /2. In the superfluid phase (|) acquires an 
expectation value. We write 
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(5) 



where % ~ ((])) and the scale m%/(27t) was introduced in order 
to have a correctly normalized boson field. In order to get 
a well defined perturbative expansion we add and subtract a 
kinetic term for the boson field to the lagrangian. We include 
the kinetic term in the free part of the lagrangian 



V2 

"V"^ /3o + 8/j + ai— + (j)o(o+ + a- 
2;7^ 



The interacting part is 

= g(^''^a+^>(p+ll.cJ+|J^>^a3'^' 

-^'('3o + |^)cp. (7) 

Note that the interacting part generates self energy corrections 
to the boson propagator Using Eq. (O we can show that to 
leading order in e these self energy corrections cancel against 
the negative of the kinetic term of the boson field in £/. We 
have also included the chemical potential term in £/. This is 
motivated by the fact that near d = 4 the system reduces to 
a non-interacting Bose gas and /j ^ 0. We will count /j as a 
quantity of (9(e). 

The Feynman rules are quite simple. The fermion and bo- 
son propagators are 



G(po,p) 
Dipo,p) 



{po + 5;u)2 - Ep 



PQ + 5/j + (ap ~x 



p2 ' 



(8) 



where &p=p'^/ (2m) and Ep= + x 

scription for the interacting theory is po 
fermion-boson vertices are igo"^ and insertions of the chemi 
cal potential are //JO3. 



The "/-delta" pre- 
*po(l+i§). The 



III. THERMODYNAMIC POTENTIAL 

The calculation of the thermodynamic potential with non- 
zero 5/j is very similar to the 8/j = calculation of Ref. 
The main difference is that the 5/j dependent pieces get con- 
tribution from momenta pin a window p < y/lmd/j iflTl [isll . 
In particular, the first one-loop diagram in Fig.[T]gives a 0(1) 
5/j contribution to the effective potential. The second diagram 
is also 0(1). However, the 5/j contribution is 0(e) since the 
diagram is proportional to an insertion of /j ^ e and there is 
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FIG. 1 : Leading order contributions to the thermodynamic potential 
in the e-expansion. The solid lines are fermions, and the dashed lines 
are bosons. The "x" represents insertion of ^703. 



no 1 /e enhancement from the finite volume loop-integral for 
momenta p < ^2m5/j. The two-loop diagram is 0(e) due to 
the factors of from the vertices. 

The first one-loop diagram from Fig. [T]gives: 



m' 



dpo d'^p log[{pQ + 5pf-El 
2% {2%Y 



(9) 



We divide the free energy by factors of m'^'+' = m^^^ to look 
at dimensionless quantities for convenience. Without the loss 
of any generality, we choose 5/j > 0. Then, 
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0(5a'-X) 
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. 2x^-(S^-X)'(5/^ + 2x)Q(5a'-X) 
24%^tii^ 



■0(e) 



with Pq = 2;7i-\/5/j2 — X^. The contribution from the second 
one-loop diagram in Fig.[T]gives 



dpo d''p Gn{po,p) -G22{po,p) 
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(11) 



0(5a'-X) 



where /j ~ 0(e). Thus the leading order effective potential is 
VQix;5p) _ a + b _ (5a(-x)^(5a' + 2x), 



24-Khn^ 
12iiV 



-0(5^ -X) 



(12) 



In Fig. 121 we plot the leading order effective potential Vo as 
a function of % for various values of 5/j. For 5/j < 2/j/e, the 
ground state is a superfluid and the potential Vo is minimized 
by xo = 2^/e. At 5^ = 2p/E, Vo{x) = Vo{xo) - mY/{3Kh^) 
for all X < S/Jc = 2/j/e. This is qualitatively different from 
weak coupling BCS results where the normal and BCS phases 
are separated by a potential barrier. At leading order in 
the e expansion the superfluid phase is stable all the way 
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FIG. 2: The leading order dimensionless effective potential V = 
3y (y^ji^ ^ function of x/fA for increasing values of 5;U. The 
dashed curve shows the result for 5^ = 0, and the solid curves are 
for non-zero 5/j. For 5^ < 2/j, superfluid phase has lower energy. 
For 8/j > 2jA, the normal phase is preferred and for b/j = l/j = 
a unique phase equilibrium characterized by the second curve from 
bottom exists. 



up to 5/j < Xq, compared to the weak coupling BCS result 
S/J ^ X0/V2. In addition to that, there is no unstable gap- 
less Sarma phase for any value of 8/j, while 
BCS calculations predict the presence of a Sarma phase for 
Xo/2 <?>iJ< xo- 

To understand the nature of the transition from superfluid 
phase to normal phase better, it is necessary to consider the 
next-to-leading order corrections. A priori it seems unlikely 
that the "flatness" of the potential for X<dfjal the critical 
will be maintained at higher orders in the expansion. The cal- 
culation can be simplified by expanding 8/j around the leading 
order value at the critical point, = 2/j/e + ^8|J'. 

The first one-loop diagram gives 



(27t)^'^ 



(13) 



r 



^ ^ 7-3yg + 31og(m7t/x) ^' 



'^+1 J i2%Y ^ e 



i^-Ep)&{2p/£-Ep) + 0{e') . 



Je ~ 0.57722 is the Euler-Mascheroni constant. The factors 
of m'^ are understood to be expanded to the appropriate order 
in e = 4 — c/. The 0(1) contribution from the last term was 
akeady calculated in Eq. ( fTOl ). The 0(e) can be calculated 
analytically but the expression is not very illuminating and we 
will not write it explicitly. 

At next-to-leading order the contribution from the second 



one-loop diagram in Fig. [T]is 



^ ^ l-2Y£ + 21og(4m7r/x) ^' 



(14) 



^\/^-X^-X^sinh-' 
0(2A'/e-x), 
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where we have used 5/i « 2/j/e + 0(e). The two-loop contri- 
bution from Fig. [T] is 



cjd+'p dd+'q Gn{p)G22{q)D{p-q) 
{2kY+^ (27i)' 
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where 



/ dx dy 
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[f{x)-x][f{y)~y] 
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is the 5^ independent piece and 
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k{x,y) + ^k{x,yf-xy 



0(5^ -X) 



0(5^ -X), 



are the 8/j corrections. We have defined the functions 

j{x,y) =f{x)+f{y) + {x + y)/2, 
k{x,y) ^f{x)-f{y)-{x + y)/2, 



(18) 



- 1 . 



Numerical evaluation gives C w 0.14424. The contribution 
from the two-loop diagram is 0(e). Therefore we can use 
5jU = 2ju/e in the integrals D, E at this order of the calculation. 

At next-to-leading order in the e-expansion, the effective 
potential is 



Vi ^i{a + b + c), 



(19) 



where we use 5/j = 2/j/e + e5/j' + (e^ ) . Results are plotted in 
Fig.[3]for various values of 8/j'. The next-to-leading order re- 
sult is qualitatively similar to weak coupling BCS theory with 
a stable and unstable superfluid phase. The stable superfluid 
phase is at Jst] 



Xo«|^[l + (3C-l+log2)e] 



(20) 
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FIG. 4: Wavefunction renormalization. Same notation as in Fig.[T] 



Fig.m The condition for finding the pole in the fermion prop- 
agator det [5^^ (/?o,p)] = now reads: 



det[G V/'o,p)-J:(/'o,p) + ;A'(J3] =0 . 



(25) 



FIG. 3: Tiie next-to-leading order dimensionless effective potential 

V' = 3^(-^) as a function of Factors of iiniifl'^ are un- 

derstood to be expanded to the appropriate order in £ = 4 — d. The 
dashed curve shows the result for 5/j = 0, and the solid curves are for 
three different values of 6/j near 2ju/e. 



The critical chemical potential 5/j£ when the normal and stable 
superfluid phase are in equilibrium can be determined analyt- 
ically by equating the corresponding pressures: 



yi(X = 0;8A(.)=^i(Xo;§A'c) 



(21) 



This is further simplified because Y\ (xo; S/Jc) = V'l (Xo, 0). We 
find 



e 



, 12C-8+51og2 ■ 
1 + ; 



(22) 



i^(l-0.4672e) 



IV. FERMION DISPERSION RELATION 

The fermion dispersion relation is determined by the pole 
in the fermion Nambu-Gor'kov propagator. From the inverse 
propagator 5"' (po,p) ~ G^^ipo,p) at leading order, we get: 



det[G-'(/7o,p)] =0, 



(23) 



which gives po = y + Xo ^ ^A*- ^.t leading order in the e- 
expansion, the quasiparticle energy has a minimum at \p\ = 
with a gap A = X() — 5iJ = 2ij/e — 5/j: 



The self-energy contribution E(/:'(vn) is diagonal in the 
Nambu-Gor'kov space and we find [511: 

Ui{p)=-8^l -^^G22{q)D{p-q), (26) 

^22{p)^-g'J j^-^Gniq)D{q-p). 

Close to the minimum of the dispersion relation, we only need 
to use the leading order values po =Xo — S/J, COp = to eval- 
uate T.{po,p) d. We write r(/?o,p) « J^'^HXo - 8a',0) + 
^*'Hxo - Sj", 0)cOp/xo- Solving Eq. (|25) gives 



po =A- 



2X0 



(27) 



y(0) , y(0) 

A=Xo-Sa^ + / "I 



coo =A' + ' 



„(0)_„(0)_„(1)_„(1) 
■^22 ^11 "^11 ^^22 



The contribution from E(/7o,p) to the dispersion relation at 
5jU = has been calculated before [5]. The 8/j dependence 
comes from momentum q integrals that involve factors of 



©(8/j — yco^+Xo)- Previously in Eq. (1211 1 we found that 
~ Xo- Thus for 5/j < S/j^, where the superfluid phase is 
thermodynamically stable, E(/?o,p) is actually 5/j indepen- 
dent. Therefore, we get 

A =^ [1 + (3C- 1 - 81og3 + 131og2)£] - 8^; (28) 



~ e 
coo =2iJ 



[l-0.345e]-8A(, 



cot 



2Xo 



(24) 



We notice that A decreases linearly with 8/j even at next-to- 
leading order. At the critical 8/ic , the gap is 



The gap decreases linearly with the asymmetry 8/j and van- 
ishes at the critical value 8/Jc = 2/j/e, where the gapless modes 
are the ones associated with the normal phase. 

At next-to-leading order, the fermion propagator S{po,p) 
gets a contribution from the self energy diagrams shown in 



6C + 2-481og3 + 731og2 _ 
Ac =p « 0.244jU , (29) 



where we have set e = 1 . Thus, at this order of the calculation 
there are no gapless superfluid modes for 8/j < 8/Jc. 
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FIG. 5: Thermodynamic potential as a function of the supercur- 
rent near the upper and lower critical chemical potential. We 
plot fl/C as a function of x for different values of h. The scal- 
ing variables x.h are defined in the text. The upper panel shows 
h = (-0.04, -0.03, -0.023, -0.01) and the lower panel shows h = 
(0.10,0.172,0.20,0.25). 



V. STABILITY OF THE HOMOGENEOUS PHASE 

We observed that there are no gapless fermion modes at 
8/j£ . Nevertheless, given the size of next-to-leading order cor- 
rections to 5/Jc and A the presence of gapless or almost gap- 
less fermions can certainly not be excluded. Moreover, gap- 
less modes might exist near the unitarity limit at finite scatter- 
ing lengths. Gapless fermion modes in BCS-type superfluids 
cause instabihties of the homogeneous superfluid phase, see 
IM m m, m m. The dispersion relation Eq. ^ is 
BCS-like, and we therefore investigate the stability of the ho- 
mogeneous phase with regard to the formation of a non-zero 
Goldstone boson current ~ V(|)/m. We note that the dis- 
persion relation is only weakly BCS-like, coq = (9(e) whereas 
Xo = 0(1). As a result any current that is formed is small, and 
we can neglect terms of higher order in v^, or inhomogeneities 
in the absolute magnitude of Xo- The effective potential for Vj 
is 



^o(v.v)=^pv.? + 1 



^E„@{-Ep) , 



(30) 



where p = nm is the mass density and Ep= A — p- + (cOp — 
a)o)^/(2xo) is the dispersion relation in the presence of a non- 
zero current. At leading order it is sufficient to compute the 
integral Eq. ( [30l ) in c/ = 4 dimensions. A rough analytical es- 
timate can be obtained by approximating the measure of the 
angular integration J^^dyy^T^-y^ 
scaling variables 



/j' I dy. Introducing the 



(31) 



_(30jt3)2p2p9^ 

(3071^)2 P^(Ao-5a<), 



with Ao = A(5;U = 0), a = con/(2mYo) and P = l/(8xo»i") we 
can write the effective potential as ll2lll 



n =yo(v,) - yo(o) = c [f,, W - /,,(o)] 



(32) 



1 



a' 



(307t3)4 p3p9 



where 



fhix) 



(h +X)V20(/J +jc) - (/l -x)V20(/i 



(33) 



The functional //, (x) has non-trivial minima in the range h S 
[—0.067,0.502]. This means that there is a range of values 
for // = — (Ao — 8/j) / Aq in which the ground state support 
a non-zero supercurrent. The size of this window is para- 
metrically small, (9(e^), and so is the magnitude of the cur- 
rent, Vs ~ e"/^. Using the leading order results for p and 
the fermion dispersion relation we get i/„,„, = —0.24 and 
Hnuix = 1 .84. Numerical results for the complete energy func- 
tional ind = 4 are shown in Fig.|5] The result is qualitatively 
very similar to the approximate solution, but the supercurrent 
window shrinks by about a factor 3. We get //„„„ ~ —0.08 and 
H„,„, = 0.63. 



VI. CONCLUSIONS 

We used an expansion around d = 4 — E spatial dimension 
to study the phase structure of a cold polarized Fermi gas at 
infinite scattering length. At next-to-leading order we find a 
single first order phase transition from a superfluid phase to a 
fully polarized normal Fermi liquid. The critical chemical po- 
tential is 5/Jc = 1 06^1. We also find an unstable gapless Sarma 
phase. We observe that <9(£) corrections are sizable, and the 
presence of gapless superfluid phase cannot be excluded. We 
show that a gapless superfluid is unstable with respect to the 
formation of a supercurrent. Recent experiments yj indicate 
the existence of at least one intermediate phase between the 
superfluid and fully polarized normal state. This suggests that 
a gapless superfluid phase or partially polarized normal phase 
is stabilized by finite temperature effects, finite size effects, or 
higher order corrections in the e expansion. 
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VII. NOTE ADDED 

Independently of this work Nishida and Son studied the 
phase diagram of a polarized Fermi liquid ifzsll . Where the 
two investigations overlap, they agree. Since these results ap- 
peared Arnold et al. computed the NNLO contribution to the 
ground state energy of an unpolarized Fermi liquid ll26ll . The 
NNLO correction turns out to be large and destroys the ap- 
parent convergence and nice agreement with Green Function 
Monte Carlo calculations observed at NLO. This is maybe 
not entirely surprising, the naive e expansion of the critical 
exponents in the Ising model also shows poor convergence 



properties at higher orders. The result of Arnold et al. im- 
plies that beyond NLO the e expansion has to be improved 
by combining the d = 4- — e expansion with information from 
the d = 2 + £ expansion, and by using Pade approximants or 
similar methods. 
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